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FOREWORD 


As  the  speed  of  air-to-water  missiles  increases,  the  water 
impact  decelerations  which  these  missiles  experience  at  water  entry 
become  increasingly  important  to  structural  design.  For  example, 
at  low  speeds  of  even  a few  hundred  feet  per  second,  decelerations 
of  several  hundred  gravitational  units  are  common.  At  these  ent^y 
speeds,  and  increasingly  more  so  as  the  speeds  Increase,  the  com- 
pressibility of  water  has  important  influence  on  the  peak  impact 
loads  and  decelerations. 

The  research  reported  here  considers  the  effect  of  water 
compressibility  in  the  fundamental  case  of  vertical  impact  of 
round-nosed  missiles,  by  applying  mathematical  techniques  pre- 
viously used  only  in  supersonic  aerodynamics. 

This  work  hcs  been  supported  jointly  by  the  Office  of  Naval 
Research  (Project  TOWER:  Treatise  on  Water-Entry  Research)  and 
this  Station  under  Local  Project  701,  from  1950  to  195^* 

This  report  is  transmitted  for  information  only.  It  does  not 
represent  the  official  views  or  final  Judgment  of  this  Station. 

It  presents  information  released  at  the  working  level  that  Is  still 
subject  to  modification  and  withdrawal. 


D.  W.  STEEL,  Acting  Head 
Research  Division 


Released  under 
the  authority  of: 

D.  W.  STEEL,  Head 
Underwater  Ordnance  Department 
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ABSTRACT 


Compressibility  effects  during  the  initial  impact  i-.iase  of 
air-to-water  missiles  are  presented  in  the  form  of  a review  of  ex- 
isting theoretical  and  experimental  work  conducted  at  the  U.  S. 
Naval  Ordnance  Test  Station  during  the  last  few  years. 

The  theoretical  section  of  the  report  presents  three  differ- 
ent approaches  to  the  problem  of  finding  the  pressures  which  act 
upon  the  missile  when  it  strikes  the  water  with  a velocity  much 
lower  than  the  speed  of  sound  in  water,  taking  into  account  the 
compressibility  of  the  water. 

The  experimental  section  of  the  report  presents  results  ob- 
tained in  measuring  the  impact  pressure  of  missiles  or  spheres 
striking  the  water  surface. 
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INTRODUCTION 


Compressibility  effects  during  the  first  instants  of  water- 
entry  impact  of  an  air-to-water  missile  are  of  interest  both  from 
a structural  and  a ballistic  viewpoint.  During  the  initial  impact, 
the  effects  of  compressibility  are  important,  whereas  during  the 
latter  phases  of  entry  the  incompressible  theory  is  adequate  for 
most  proDlems.  It  should  be  emphasized  that  the  amount  of  experi- 
mental data  concerning  the  compressible  phenomenon  is  very  limited, 
so  that  the  treatment  of  compressible  effects  is  essentially  a 
theoretical  one  at  the  present  time. 

Extremely  high  pressures  can  be  developed  on  the  nose  of  a 
missile  at  water  impact.  If  a flat,  rigid  disk  is  considered  as 
striking  a water  surface  normally,  the  instantaneous  pressure  due 
to  the  compression  wave  set  up  in  the  water  is  p »£>cV  (Appendix  A) 
v are  p is  the  density  of  the  water,  c is  the  velocity  of  sound  in 
* <s  water,  and  V is  the  normal  velocity  of  the  disk.  If  this  prob- 
lem is  considered  from  the  incompressible  viewpoint,  the  pressure 
ao  impact  has  an  infinite  value.  In  the  case  of  a sphere  striking 
a water  surface  normally,  there  is  a small  area  over  which  such  a 
shvck  will  exist  for  a finite  interval  during  which  the  intersec- 
tion of  the  sphere  with  the  water  surface  expands  faster  than  the 
velocity  of  sound  in  water.  For  low  entry  velocities,  V<JCc,  the 
time  during  which  the  rate  of  growth  of  the  wetted-surface  radius 
is  supersonic  is  very  small.  Hence  for  most  cases  of  water  entry, 
the  impact  or  shock  phase  occurs  during  a very  short  interval, 
usually  of  the  order  of  a few  microseconds.  This  factor  alone 
makes  experimental  work  very  difficult  to  perform  since  instruments 
with  very  rapid  response  are  required  in  order  to  measure  events 
occurring  during  a short  interval.  In  addition,  there  are  indica- 
tions of  air  trapped  between  the  missile  and  the  water  which  acts 
as  a cushion  in  its  effects  on  the  shock  phase.  Since  the  pres- 
sure plugs  in  the  missile  have  a finite  diameter,  any  pressure 
measurement  indicates  an  averaged  value.  The  phenomenon  of  trapped 
air  may  also  delay  and  reduce  the  peak  pressures. 

The  structural  effect  of  the  shock  phase  is  of  interest  in 
connection  with  the  deformation  of  the  nose  of  the  missile,  but 
it  is  not  well  known.  The  structural  effect  is  also  of  theoretical 
interest  since  it  combines  hydrodynamic  and  elasticity  theory  ir 
the  solution  of  the  problem. 
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Ac  -ono  part  oi  tin  over-all  problem  of  water  entry  shock,  the 
Impact  phase  In  important  since  it  is  the  first  to  occur  and  the 
subsequent  behavior  of  the  missile  is  dependent  on  what  happens 
during  this  phase.  The  fact  that  the  compressibility  effects  are 
more  pronounced  during  this  phase  means  that  only  theoretical  and 
experimental  work  of  a fundamental  nature  can  revaal  their  relative 
importance. 

The  incompressible  theory  of  water  entry  of  blunt  bodies  has 
been  extensively  developed  by  Shiffman  and  Spencer  (Ref.  l)  and  by 
others.  Based  on  this  theory,  the  pressure  distribution  on  a 
sphere  entering  water  vertically  has  been  determined  at  the  Naval 
Ordnance  Test  Station  (NOTS),  and  the  theoretical  results  have  been 
checked  experimentally  through  pressure  measurements  with  piezoelec- 
tric gages  at  the  stagnation  point  of  a 12-inch  ’diameter  sphere. 

Incompressible  theory,  however,  predicts  an  infinite  pressure 
at  the  stagnation  point  of  a blunt  body  when  it  comes  in  contact 
with  the  water  surface.  Actually,  it  is  expected  that  during  the 
impact  phase,  pressures  of  the  order  of  the  shock  or  piston- 
impulsive  pressure  (?cV  are  experienced  for  a few  microseconds. 

Since  a theoretical  prediction  of  the  pressures  during  this  phase 
requires  that  account  be  taken  of  the  compressibility  of  the  water, 
the  formulation  of  this  problem  and  several  solutions  will  be  con- 
sidered. The  problem  was  first  stated  by  L.  Trilling  of  this  Sta- 
tion, who  based  his  formulation  on  the  theory  of  weak  waves  and 
obtained  a solution  to  the  two-space  dimensional  problem  after  mak- 
ing several  assumptions.  The  problem  was  later  taken  up  at  NOTS  by 
R.  H.  Owens,  who  reformulated  the  problem  using  the  more  general 
Reisz  method,  but  maae  no  attempt  at  a detailed  solution.  The  ver- 
tical .mpact  of  a sphere  on  the  water  surface  was  considered  at 
NOTS  by  R.  H.  Korkegi,  using  the  retarded  potential  solution  as 
suggested  by  P.  E.  Marble,  consultant  to  this  Station,  from  which 
actual  pressure  distributions  were  obtained. 


TRILLING'S  SOLUTION  FOR  TWO-SPACE  DIMENSIONS 


Since  this  report  deals  with  the  first  contact  of  a striking 
body  with  a plane  water  surface  where  the  compressibility  of  the 
uater  is  to  be  considered,  the  problem  resolves  itself  into  one 
concerning  the  propagation  of  a pressure  wave  of  finite  amplitude 
into  the  water.  This  involves  considerations  in  terms  of  the  dy- 
namics of  a compressible  liquid. 

_ When  a body  strikes  a water  surface  at  a velocity  V where 
|VJ/6«1  (c  is  the  velocity  of  sound  in  water),  the  theory  of  weak 
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waves  (Ref.  2)  may  be  applied  to  compute  its  motion.  The  flow  is 
Irrofcaticnal  and  isentrcpic,  and  the  velocity  potential  satisfies 
the  wave  equation  up  to  terms  of  the  order  of  V2/c2.  The  boundary 
conditions  are  applied  on  the  undisturbed  free  surface,  and  the 
displacement  of  the  body  is  neglected.  The  effect  of  splash  is 
ignored. 

The  velocity  potential  f satisfies  the  equation 

(1)  Vtl  = o 

The  velocity  V and  the  pressure  p are  determined  by 

(2)  v = vv  , p = -e^t 

Just  before  impact,  the  fluid  is  at  rest  or 

0)  <P(  o,x>y>-i)^o 

Since  the  differential  equation  and  the  boundary  conditions  have 
been  linearized,  the  solutions  for  the  vertical  and  horizontal 
components  of  motion  may  be  found  separately.  In  this  report  only 
the  vertical  component  will  be  considered. 

The  boundary  condition  for  the  vertical  component  is 

(fc)  P( tj  X,0,1)  - 0 o-»  S' 

a)  = v ov>  5 

where  S is  the  wetted  surface  of  the  body  and  S'  is  the  horizontal 
free  surface  of  the  liquid.  The  problem  (Eq.  1 - 4)  is  similar  to 
that  of  a lifting  three-dimensional  wing  in  a steady  supersonic 
stream  for  the  case  of  two-space  dimensions  (Ref.  3>  p.  73).*  and 
the  methods  of  supersonic  airfoil  theory  may  be  adapted  to  solve  it. 

When  a two-dimensional  body  strikes  a plane  water  surface  at 
a constant  velocity  V,  its  wetted  surface  jx/<xe(t),  grows  at  a 
rat'  proportional  to  the  slope  of  the  body  section  since 
dxe/dy  = xe/v.  The  potential  function  satisfies  the  equation 

• (5)  (fu-  ^ ( Wcx  f fty) 

with  the  initial  condition 

= o 
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The  vertical  component  of  motion  satisfies  the  boundary  conditions 


K > Xe(t) 


Py(*,X,0)=V  > 


x < Kelt) 


Figure  1 shows  the  growth  of  the  strip  xe(t)  in  the  plane  y 


= 0. 


r. 


x 


The  lines  AB,  AB1  are  the  traces  of  the  characteristic  Mach  cone 
tangent  to  x(t)  in  the  plane  y = 0.  Since  the  equation  of  motion 
is  hyperbolic  with  constant  leading  coefficients,  any  disturbance 
travels  in  the  fluid  field  at  the  constant  velocity  c.  In  the  re- 
gion B'  OB,  the  edge  of  S moves  outward  at  a velocity  xe>c,  and 
therefore,  in  the  region  BOB'C,  no  point  on  the  body  surface  can 
be  influenced  by  the  fact  that  the  body  has  finite  width.  The 
pressure  in  that  region  can  be  determined 
wave  theory  (Ref.  3>  P*  4l)  and  is 

and  since  V<<(c  and  xe>c  for  this  domain, 
proximated  very  closely  by 

poc  (O  c V 

To  investigate  the  flow  in  regions  OBC,  OB'C^it  is  convenient 

to  approximate  the  curve  OB  by  a straight  line  whose  angle  & with 

the  horizontal  x-axis  is  small,  since  tan^i  ^ -3L.  The  cuive  BD  is 

1 c 


from  one -dimensional 

; V 

the  pressure  can  bo  ap- 
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approx lmu I O'.'  by  a straight  line  whose  angle  with  the  vertical  (ct) 
ex's  la  0 ..  The  point  P Is  any  convenient  point  on  x0(t).  The 
origin  of  the  coordinate  systom  io  now  shifted  to  B.  The  flow  In 
if-i-  r • .-.ion  CBC  is  conical  in  the  sonoo  Introduced  by  Buoomanri 
i(  >i  ’i).  It  Is  char .cterlaod  by  th«*  absonoc  of  any  length  parun*- 
iti  in  the  o^uatlon  of  motion  «■  the  geometry  of  the  field-  The 
- i •»*  Uy  and  pressure  depena  only  on  tin*  parameters  £ = x/ct  and 
. ’ In  the  f plane,  the  Much  h.m:o  from  B appourn  ao  the 

i*  iroln.  The  boundary  BD  io  a segment  ■,»  th"  f uxln  The 
*„  ltf  contact  of  the  Mach  cone  from  B with  the  plain*  wuvoa  from 
the  points  PP1  or>  t;*«  unit  circle  {Pig  ;*),  s 1 1 ua « e<l  ut 


Cos  y7 » 


VI, 


Tho  differential  equation  f«r  *h  pressure,  >i  t»in  I in 


troduclng  p(  ^ ,b)  Into  tho  wavo  oqual  *n,  jgq. 
to  polar  coordinates  vr , Y In  tho  ] fj  plane  is 


;i*  i * i i»  - I ■ 


(i-‘r‘)P^t-  ~~  fir  ' 'a  Pr, 
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The  angular  coordinate  <T" is  now  left  unchanged,  while  the  radial 
coordinate  is  transformed  as  suggested  by  Tchaplygine  (Ref.  J>, 

p.  28)  

/ — / / —cr& 

(6)  5 * 


The  unit  circle  and  the  origin  remain  invariant  and  rays  from  the 
origin  are  transformed  into  themselves,  although  not  point  by  point. 
In  particular,  the  edge  D which  was  situated  at  the  pointtTD  = tan<%> 
1 * n.iw  defined  by 

Sp  = c.0t  4 (/-, 

In  the  S,  tin,  pressure  and  the  velocity  components  satisfy 

file  equation 

/\,  <-  A f ,0 

The  filin’ I lens  |>,  u,  v are  t.horofore  real  parts  of  analytic  functions 
1* , 11,  V of  Mi*  complex  variable  Se^®.  Since  v is  constant 
along  vnnl.ihep  there  It  follows  from  the  momentum  equation 

tlm!  <>.  but  a ong  the  rooJ  axis 

d/'  - -L  $£- 

o»  y 0~  o>  S’ 

uml  because  of  the  invariance  of  angles  undor  the  Tchaplygine  trans- 
formation^? vanishes.  From  the  Cauchy-Riemann  conditions,  it 
follows  tha*  the  harmonic  conjugate  p*  of  p is  constant  along  CD 
>i ml  may  be  set  equal  tc  zero.  On  F1  CF  where  / y/<"  f , no  disturbance 
<..ii  affect  the  pressure  field,  so  that  the  pressure  vanishes  there. 
Since  the  pressure  is  continuous  in  the  field,  the  pressure  along 
FCF'  the  same  as  behind  the  wayes  from  OB,  that  is  p =£>cV  along 
’’F  unit  -£>cV  along  CF'.  To  simplify  conditions  along  the  real  axis, 

It  In  convenient  to  introduce  the  homographic  transformation 

in  vV  — 

/-  So 


which  leaves  the  unit  circle  and  the  real  axis  invariant,  and  puts 
the  point  D at  the  origin.  The  points  FF1  are  now  located  at 

= £ Q-n  L -j/  /—  ■=.  EL~~  ^ 

Along  the  positive  real  w-axis  p is  zero,  and  since  it  is  harmonic 
it  must  be  antisymmetric  with  respect  to  the  real  axis.  With  re- 
spect to  the  negative  cal  w-axis,  it  is  symmetric  by  virtue  of 
Schwartz's  reflection  principle,  since  p*  vanishes  there.  The 
w-plane  (Fig.  3)  is  therefore  a double-sheeted  Riemann  surface  with 
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a cut  along  CD.  A reflection  of  P(w)  in  CD  defines  P(w)  in  the 
lower  sheet.  The  solution  in  the  physical  plane  has  a discontinu- 
ity along  CD,  where  the  body  can  support  a pressure  discontinuity, 
but  it  is  continuous  in  the  fluid,  along  the  positive  real  axis. 


If  the  douole-sheeted  surface  is  unwound  by  the  transformation  t 
the  points  F,  F'  are  mapped  into  four  points  \ & £ gp 0 'w , 
and  the  boundary  conditions  for  P $)  are  given  on  the  unit  - \ 
circle  as  follows 

!UPC1)=0> 

Re  Ptt)^ecV  j €£p5? 

W Cl 

R&  pcx)*-e*v  j < *ra>  < 


The  function  which  satisfies  these  conditions  is  (Ref.  3,  p. 


pc»  = 


r ^ Ci-WO-z,) 


79) 


The  pressure  distribution  on  the  body  surface  is  found  by  retracing 
back  through  the  transformations  in  Eq.  8,  7,  and  6.  The  result 
for  small  $ (tan  9^9,  ) is 


P * [A  0 i-  B &,)] 
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where 

* ' /-S  1 Z ()+$%)( H-Sf) 


) 


Sp  — cot  &z  (/- y /-  ta-j'6^ ' ) } 


cr  - 


J£_ 

Or 


s-  tJE^s 

cr  ) 


Several  factors  restrict  the  usefulness  of  this  solution. 
First,  che  boundary  of  the  body  has  been  approximated  by  straight 
line  segments.  For  this  approximation  to  be  reasonable  it  is  nec- 
essary to  make  these  segments  very  short.  In  this  connection 
Evvard's  method  (Ref.  5)  based  on  the  appropriate  fundamental  solu- 
tion of  the  wave  equation  might  be  used  in  order  to  wliminate  ap- 
proximation of  the  boundary  by  straight  line  segments.  Coupled  to 
this  restriction  is  the  fact  that  to  extend  the  solution  any  fur- 
ther requires  the  use  of  complicated  methods  of  superimposing  other 
solutions  and  then  matching  boundary  conditions.  Because  of  the 
complexities  of  such  a solution  and  the  questionable  accuracy  of 
the  solution  already  at  hand,  such  an  attempt  does  not  seem  to  be 
justified.  Since  the  two-space  dimensional  case  is  essentially  of 
academic  interest,  it  is  worth  noting  that  the  problem  can  be 
treated  by  using  the  fundamental  solution  of  the  wave  equation. 

This  has  been  done  by  R.  H.  Owens  using  Rle3z's  method  (Ref.  6). 

It  should  be  noted  here  that  the  retarded  potential  cannot  be  for- 
mulated in  two  dimensions  because  of  Huygen's  principle. 

Riesz's  method  is  presented  in  Appendix  B.  The  advantage  of 
the  method  is  that  it  clearly  shows  what  boundary  and  initial  con- 
ditions are  necessary  to  solve  a problem  in  hyperbolic  differential 
equations.  Recalling  that  in  m dimensions  the  hyperbolic  distance 
R is  gi  by 

& *-  Z (t-  *? 

j- 1 

and  since 

/0,  = o £J  — P A/e'Yriber'ti&n'^ 

R2-m  satisfies  this  m-dimensional  wave  equation.  However,  an  at- 
tempt to  build  up  a solution  irom  this  so-called  fundamental  solu- 
tion fails  when  m >3  because  any  integral  involving  R2"1”  diverges. 
However,  an  integral  of  RO.+2-ro  where  a>  m - 4 will  converge  (the 
use  of  analytic  continuation  is  the  basis  for  Riesz's  method)  and 


8 
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its  limit  as  a^O  may  be  considered  as  the  analytic  continuation  of 
the  integral  of  Ra+2-m  into  the  domain  0^.a^m  - 4,  i.e.,  it  de- 
fines the  integral  of  R2_m.  In  order  to  perform  this  analytic  con- 
tinuation it  ib  necessary  to  integrate  by  parts.  However,  the  in- 
tegral is  multiple,  and  Sreen's  theorem  is  used,  which  corresponds 
to  integration  by  parts.  For  m = 4 the  familiar  retarded  potential 
solution  is  obtained.  However  this  can  be  done  more  simply,  and 
constitutes  the  basis  for  the  formulation  of  the  problem  by  R.  H. 
Korkegi . 


KCRKEGI'S  SOLUTION  FOR  THREE-SPACE  DIMENSIONS 


The  formulation  of  the  problem  is  again  essentially  that  of 
L.  Trilling.  The  first  few  instants  of -water  entry  during  which 
only  a very  small  part  of  the  sphere  is  submerged  will  be  consid- 
ered. The  depth  of  penetration  is  small  compared  to  the  radius  of 
the  wetted  area;  hence,  to  the  order  of  approximation  of  this  anal- 
ysis, boundary  conditions  will  be  satisfied  in  the  plane  of  the 
undisturbed  water  surface  rather  than  on  the  curved  surface  of  the 
sphere.  Since  the  surface  of  contact  of  sphere  and  water  is  the 
only  disturbance  present,  the  analytic  problem  is  that  of  deter- 
mining the  flow  field  that  is  due  to  an  expanding  disk  of  disturb- 
ances. In  addition,  the  problem  will  be  restricted  to  considering 
tile  time  during  which  the  rate  ,f  expansion  of  the  wetted  surface 
exceeds  the  rate  of  wave  propagation  in  water  and  the  free  surface 
of  the  water  beyond  the  area  of  contact  is  undisturbed.  The  fol- 
lowing assumptions  are  made: 

1.  That  the  flow  Is  irrotational  and  isentropic. 

2.  Th»t  the  velocity  of  the  body  V and  the  rate  of  wave  prop- 
agation in  water  c are  constant. 

3.  That  V/c  is  much  smaller  than  unity. 

4.  That  boundary  conditions  are  satisfied  in  the  plane  of  the 
undisturbed  water  surface. 

5.  That  the  sphere  is  considered  rigid. 

The  condition  of  irrotationality, VxU  = 0,  is  identically  satisfied 
by  » potential  such  that  n r=7j^.  (U  is  the  velocity  vector  in  the 
flow  field).  With  a system  of  coordinates  x,  y,  rt  (or  r,  a,  It) 
fixed  in  space,  the  analytic  problem  consists  of  vhe  following  par- 
tial differential  equation,  boundary*  and  initial  conditions.  Fig- 
ure 4 shows  the  area  of  contact  of  sphere  and  water  represented  by 
a disk. 
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(8) 

V2<f-  l/c z 

(9a) 

if(  r,o,  b)  -€>  > 

r > Zt'b) 

(9b) 

b)~£>  ) 

r > Hit) 

(9c) 

r < R.lt ) 

where 

n l,p)  = (r} 

n 

“N 

f u 

rz  t 2*  ^0 


R(t)  = the  instantaneous  radius  of  the  disk 
^J/,2-  = coordinates  of  the  field  point 
f,  = coordinates  of  a variable  point  in  space 
t~  = fixed  time 
"2^  = variable  time 


PIG.  4. 


A fundamental  solution  of  the  wave  equation  (Eq.  8)  is  the 
simple  source 

¥{*},*)**  - 


n 
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It  represents  the  instantaneous  potential  at  time  t at  a 
'■’i'-’ld  point  a distance  r,  away  from  a source  of  strength  f.  Since 
the  wave  equation  is  linear,  solutions  can  be  constructed  by  super- 
imposing simple  sources.  In  particular,  the  solution  for  a dis- 
tribution of  point  sources  in  the  plane = 0 is 


V(x,y,i,b)  ~ fj fSIdkfLL. 

4 ' 


where 


It  is  found  that  the  instantaneous  pulse  strength  f is  directly 
proportional  to  the  instantaneous  and  local  value  of  This  is 
shown  by  Lagers trom  (Appendix  C)  and  yields 


H 1, 1,  °>  *)  = '£(t  1>  a>  x) 


Thus,  the  solution  can  be  written  in  terms  of  a retarded  potential 
as  follows 


(10a) 


or 


(10b) 


where 


-±  rr 

£ry' 


ct  zjt 


fid  i)Mr, 


The  limits  of  integration  of  Eq.  10b  are  obtained  as  follows: 
if  there  are  no  disturbances  prior  to  time  X = 0,  then,  at  time 
X - t,  the  field  point  P (Fig.  5),  is  influenced  by  disturbances 
occurring  within  a sphere  of  radius  ct  about  P;  since  disturbances 
are  limited  to  the  plane  t\  - 0,  the  region  of  integration  is  the 
circular  area  defined  by  the  intersection  of  the  sphere  ct  with 
the  plane  ^ = 0.  The  integrand  of  Eq.  10a  and  b is  determined  by 
the  boundary  conditions  of  the  problem. 
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2-4 

FIG.  5- 


For  a sphere  of  radius  a entering  water  with  a vertical  veloc- 
ity V the  inrtar.taneous  radius  of  the  wetted  surface  for  small 
depths  of  ptnccration  is  approximately  (see  Fig.  6) 


The  rate  of  growth  of  the  radius  of  the  wetted  area  is  hence 


When  the  rate,  of  growth  becomes  sonic  R(tc)  = c,  the  time  at  which 
the  rate  of  growth  is  sonic  is 

L _ ay 
(13)  !'c  " 
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■ >SS2a 


Since  this  analysis  will  hold  for'  times  during  which  the  rate  of 
growth  is  supersonic,  it  must  be  required  that 


t<tc  , K/t)£  fl(k)  = 


It  is  convenient  to  use  dimensionless  variables  (denoted  by  primes) 
as  follows 


(14a)  RUt)  = - I 

(14b)  t}  « £-/|^  - i 


The  dimensionless  rate  of  growth  becomes  the  Mach  number 


2 


Now  the  boundary  conditions  (Eq.  9)  are  applied  to  the  re- 
tarded potential  solution  (Eq.  10a)  in  order  to  determine  the  po- 
tential at  an  arbitrary  field  point.  Since  the  problem  is  axially 
symmetric,  no  generality  is  ’"'st  by  choosing  a point  in  the  verti- 
cal plane  y = 0.  Because  of  ;he  time  dependency  involved,  not  all 
pulses  emitted  in  the  circular  region  of  integration  illustrated 
in  Pig.  5 are  felt  at  the  field  point  P at  time  t.  Hence  it  is 
necessary  to  determine  the  region  of  integration  within  the  cir- 
cular area  of  radius  i/  c2^-  2 , with  pulses  of  strength  = V 
(Eq.  9c)  which  will  contribute  to  the  potential  at  P (x,  0,  b,  t). 
This  region  is  bounded  by  the  locus  of  the  points  of  Intersection 
of  the  instantaneous  Mach  cone  (see  Fig.  Ja)  with  the  expanding 
disk  (representing  the  wetted  surface  of  the  sphere)  for  all  times 
between^  = 0 and  ¥ = t - b/c.  Prom  Pig.  7a  it  can  be  seen  that 
pulses  emitted  for  0 <¥<£!'/,  and  ¥/>2zhave  not  yet  reached  P at 
time  t;  only  those  emitted  for  % < Z<  2ghave  influence.  Hence,  the 
area  in  the  plane  z = 0 in  which  disturbances  emitted  will  reach 
the  field  point  P(x,  0,  s)  at  time  t is  bounded  by  the  curve^(^) 
(Pig.  7b)  given  by  the  equation  1 

RHt)~  V+f 

where  

H2(Z)~  g)  = ) 


13 


:<  ■rrrr-’ViZ 


NAVORD  REPORT  3523 


This  equation  can  be  expressed  in  dimensionless  form  as  fol- 


lows 


t ]/ (X  - ? J)z  2,Z  = / X -f~ 


where  t 


,=^/4y q an<*  sPace  Quantities  are  reduced  to  dimension- 
less form  through  division  by  , that  is, 

Solving  for  % (£)  the  equation  is  obtained  for  the  curve  bounding 
the  area  of  integration  as  follows 

(15!  l'(i‘U 


where 


f-.(V)--  2 /■  i1'-  r2 

i'y] 

Hence,  from  the  retarded  potential  solution  of  the  wave  equation 
(Eq.  8)  and  the  boundary  conditions  (Eq.  9b  and  c),  the  potential 
at  an  arbitrary  field  point  (x,0,e)  at  time  t,  given  in  terms  of 
dimensionless  coordinates,  is  , ,/  / 

f*  J w 

( 1 6 ) m,  ^ 2,  t)  - -==~X^ 


S’ 


A,/  *'  <?  / J>/ 

where  £ <^)  is  given  in  Eq.  15  and  the  limits  7,  and  3,  are, the 
real  roots  of  the  quartic  equation  obtained  by  setting  fg  (f)  = 0. 
In  dimensional  coordinates  the  potential  is 

r r W J, 

(17)  7 f = 

rl  y*  i(*-  V 

where 


/,/ 


?/()=  **Vt-f*-  a 


The  region  of  influence  of  disturbances  due  to  contact  of  the 
sphere  with  the  water  surface  is  sketched  in  Pig.  8.  The  dashed 
line  represents  the  front  of  the  compression  wave  moving  out  into 
the  fluid  from  the  disturbed  region. 


14 


NAVORD  REPORT  3523 


\ 


\ From  Eq.  17 


Wave  Front 


\ Region  of  influence  ' 


/ 


/ 


0 


\ 

FIG.  8. 


*J  ? 


Due  to  the  complexity  of  its  limits,  the  double  integral  for 
the  potential  (Eq.  17)  does  not  lend  itself  to  an  exact  solution 
for  an  arbitrary  field  point.  For  x = 0,  however,  the  potential 
and  the  pressure  can  be  readily  evaluated.  In  this  case  the  area 
of  integration  degenerates  to  a circle  about  the  origin.  The  result 
is  given  in  terms  of  dimensionless  space  and  time  coordinates 

(18) 

Fortunately,  the  integral  for  the  potential  can  be  simplified 
for  the  evaluation  of  the  pressure  distribution  on  the  wetted  sur- 
face of  the  sphere.  Since  the  integrand  is  a continuous  function 
of  f and  ^ for  every  value  of  z,  the  potential  in  the  plane  z 0 
can  be  directly  evaluated. 

Now  polar  coordinates  are  introduced  in  the  xy-plane  (see 
Pig-  5)  with  the  field  point  P(r,0, t)  as  origin  . 

C = yr(x:f)2  ) e = (t^7/ 


Since  r,  — > (P  as  e — > 0,  the  retarded  potential  solution 
(Eq.  10b)  becomes 

c£  air 


(19) 

where 


Wn  zfj*  J % (<°j  e/L)J*cJe 


r z~  x 2 f-  •h-‘ 
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Because  of  axial  symmetry,  the  potential  is  independent  of  a (see 
Pig.  4).  The  domain  of  dependence  of  P(r,0,t)  in  the  xy-plane  is 
the  area  bounded  by  the  circle  of  radius  ct  about  P.  However,  as 
for  the  case  of  an  arbitrary  field  point,  only  those  pulses  inside 
the  area  defined  by  (5)  will  contribute  to  the  potential  at 
p(r,0,t).  (Figure  9 shows  the  area  of  integration  for  the  point 
P(r,0)  at  time  t.) 


y 


FIG.  9- 

This  can  be  expressed  mathematically  as  follows 

(a.)  %)=  v ,o<e^e,(e>) 

<e*(e,  e,t  - e/c)  = <<?<  ct 

Hence,  Eq.  19  becomes 

2/T 

(21)  <f(r,  Oj  t)  = -X-  f e,(B)Xe 


Jp 
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The  equation  giving  the  boundary  of  the  area  of  integration  is 
(see  Pig.  9) 

R.\t)  = siv>£&  -h  (r  -het(e)ci>se)z 
= €/&)  a ar fye>)cDsa  j-r2- 

Since  the  square  of  the  instantaneous  radius  of  the  disk  (wetted 
surface)  is 

ft*cz)  - ? 'Z'V'b  - (t- 

In  dimensionless  coordinates  is  obtained 


(22)  *C&)  ■=.-(!+■  r'cesfr)  -hjj/i-r'ces&f-h  (t-K  z 

where 

/ A/sJ'  . r*-  > 

Since  for  supersonic  expansion  of  the  disk,  (?t  for 

0~&-£TT>  i.e.,  the  area  of  integration  always  encloses  the  point 
P(r,0,t)  on  the  disk. 

For  the  two  limiting  cases  of  a point  first  at  the  center  of 
the  disk  (r1  = 0),  and  second  on  the  edge  (r1  = t,1),  the  area  of 
integration  degenerates  to  a circle  of  radius  f = -1  + y 1 + t1 ' 
and  then  vanishes,  respectively.  f 

The  equation  for  the  potential  at  a point  r on  the  disk,  given 
in  terms  of  dimensionless  quantities  is  hence 

ar 

(23)  Wr^t)  = ^ ^/-^y [(Itr'cesefi-  (*'- 

To  the  order  of  approximation  of  this  analysis,  the  pressure  is 
given  by 


(24) 


p = e % = - e 
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Hence,  the  pressure  corresponding  to  the  potential  of  Eq.  23,  and 
given  in  terms  of  the  shock  pressure,  is 


PCrjQ.t)  !_  f 

V IT  JDij - -h(t'-r,z) 


This  is  a complete  elliptic  integral  of  the  first  kind,  which  can 
be  reduced  to  the  Legendre  standard  form  (see  Appendix  D)  yielding 


(26) 


e cv 


- A(r'b')  /([?c(r}t:')J 


where 




IT [(It  t')’--  (ar’f] 


/ -h  t,~Zr,z  7)^ 


Values  of  the  elliptic  integral  )((a)  are  tabulated  by  Jahr.ke  and 
Emde  (Ref.  7). 

The  time  history  of  the  pressure  distribution  on  the  sphere 
during  impact  is  plotted  in  terms  of  the  dimensionless  radius  of 

the  wetted  area  r1  = (actually  the  projection  of  the  wetted 

c 

area  in  the  plane  of  ttie  undisturbed  water  surface),  with  the  di- 
mensionless time  t '-V<xV  as  a parameter  in  Pig.  1. 

ZCz 

It  is  to  be  recalled  that  this  analysis  is  valid  only  for 
times  during  which  the  rate  of  expansion  of  the  wetted  area  of  the 
sphere  exceeds  the  speed  of  wave  propagation  in  water,  i.e.,  for 

t < -^tj  or  t'  < 1.  For  this  case  the  free  surface  of  the  water 
2c^ 

beyond  the  area  of  contact  is  undisturbed.  If  the  expansion  rate 
were  less  than  the  speed  of  wave  propagation,  compression  waves 
would  propagate  outward  from  the  edge  of  the  area  of  contact  and 
disturb  the  free  surface  inside  the  wave  front  (mathematically 
this  means  that  ^ ^ P for  R(t)<  r<|^  + ct). 
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From  incompressible  theory  the  pressure  at  the  stagnation 
point  of  a sphere  is  given  in  terms  of  the  dynamic  pressure  by 
(Ref.  1) 


As  a function  of  the  dimensionless  time  parameter  t'=  the 

incompressible  pressure  can  be  expressed  in  terms  of  the  shock 
pressure  as  follows 


(27)  — t— 

e C V 

From  Eq.  25  the  compressible  stagnation  pressure  (r1  = 0)  is 


V ■/  !+tp 

Equation  27  clearly  shows  the  nature  of  the  singularity  of  the 
Incompressible  pressure  foraula  at  the  instant  of  impact  (f  = 0) 
while  Eq.  28  indicates  that  the  compressible  pressure  formula  has 
the  finite  value  pc/. 

Equations  28  and  29  are  plotted  in  Fig.  11,  with  a predicted 
curve  for  the  compressible  pressure  when  the  rate  of  expansion  of 
the  wetted  surface  of  the  sphere  is  subsonic.  With  increasing 
time  or  a decreasing  rate  of  expansion,  it  is  expected  that  the 
compressible  and  Incompressible  pressures  approach  each  other. 

When  the  rate  of  expansion  of  the  wetted  surface  becomes  much 
smaller  than  the  speed  of  wave  propagation  in  water,  the  effects 
of  compressibility  are  negligible,  hence,  incompressibility  theory 
is  quite  adequate  for  pressure  predictions. 


EXPERIMENTAL  WORK 


The  theoretical  work  at  this  Station  has  been  concerned  with 
the  entry  of  a rigid  body  into  water  and  the  subsequent  motion  of 
the  fluid,  taking  into  account  the  compressibility  of  the  water. 

No  account  of  the  fact  that  an  actual  object  is  not  strictly  rigid 
but  deformable  has  teen  taken  In  any  of  the  theoretical  work.  Yet 
the  structural  effects,  such  as  deformation  of  the  nose  at  impact 
or  even  failure  at  high  entry  velocities,  is  of  extreme  practical 
importance. 
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0 1 2 t,  3 4 5 


PIG.  11.  Comparison  of  the  Compressible  and  Incompressible 
Pressure;  Time  History  at  the  Stagnation  Point  of  a Sphere. 
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The  difficulties  which  arise  in  any  experimental  efforts  to 
isolate  the  compressibility  effects  are  fairly  well  recognized. 

The  fact  that  the  missile  acts  as  a compressible  medium  has  been 
noted  above.  It  is  also  true  that  the  water  surface  is  not  a 
strictly  smooth  plane.  In  addition,  there  is  evidence  that  air  is 
entrapped  between  the  bo  y and  the  water.  All  these  factors  make 
it  difficult  to  separate  the  true  compressibility  effect  of  the 
water  and  modify  all  experimental  results.  Most  of  the  experimen- 
tal work  yields  order-of -magnitude  results  which  may  act  as  guides 
to  the  actual  design  of  missiles  and  future  theoretical  ana  exper- 
imental work. 

The  normal  component  of  force  at  the  point  of  impact  for  a 
blunt-nosed  missile  entering  the  water  at  an  angle  9 arises  from 
the  pressure f'cV  sin^  where  V sin  9 is  the  velocity  of  the  missile 
normal  to  the  water  surface.  In  order  to  study  the  peak  impact 
pressures  when  the  nose  of  a projectile  first  comes  in  contact 
with  the  water,  a set  of  pressure  plugs  were  distributed  over  the 
nose  of  a hemispherically  shaped  missile  (Ref»  8).  The  plugs  arc 
thin  phosphor-bronze  diaphragms  supported  peripherally  on  an  accu- 
rately reamed  shoulder.  For  entry  velocities  cf  500  fps  or  more 
the  plugs  used  were  0.02  inch  thick  and  0.250  inch  in  diameter. 
Application  of  a sufficiently  large  pressure  gives  the  diaphragm 
a permanent  set.  The  pressure  plugs  were  calibrated  statically. 
Figure  11  of  Ref . 8 shows  a plot  of  the  pressure  at  the  impact 
point  as  a function  of  velocity.  The  agreement  with  that  calcu- 
lated by  the  formula pcV  sin  9 is  very  good,  considering  the  fact 
that  the  static  calibration  of  the  plugs  was  used-  It  Is  true 
that  the  natural  frequency  of  the  membrane  is  theoretically  high. 
For  these  membranes  the  frequency  should  be  25,000  cps.  The  eff3C- 
tive  frequency  as  an  inelastic,  deformable  membrane  was  not  meas- 
ured, so  the  true  dynamic  response  is  not  known  but  is  estimated 
to  be  about  5,000  cps.  For  a hemispherical  head  with  a radius 
a of  1 foot,  entry  velocity  V of  500  fps,  and  c of  4,800  fps,  the 
duration  of  thepeV  pressure  is  roughly  twice  the  time  at  which 
the  growth  of  the  radius  of  the  wetted  area  of  the  sphere  is  just 
sonic  or  t = iZ,  =22  times  10“®  second.  This  is  because  of  the 

reflection  oi  the  pressure  wave  from  the  point  where  the  velocity 
of  the  contact  point  is  just  equal  to  c.  The  analogy  for  the 
pressure  plug  is  a spring-mass  system  which  can  move  only  in  one 
direction  just  as  though  the  mass  were  restricted  by  a rachet. 

Hence  the  time  for  one-fourth  cycle  of  the  motion  of  the  membrane 
needs  to  be  considered.  For  5>000  cps  this  time  is  t = 50  micro- 
seconds, as  compared  to  a duration  ofpcV  of  t = 22  microseconds. 
Therefore  the  plug  should  not  be  deformed  as  much  as  the  static 
calibration  would  indicate.  However,  as  an  order  of  magnitude, 
the  check  is  very  good. 
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A laboratory  setup  for  obtaining  pressure-time  measurements 
has  been  established  at  NOTS  by  C.  R.  Nisewanger  (Ref.  9)* 
Nisewanger's  method  consists  of  making  direct  pressure  measure- 
ments on  bodies  by  moans  of  electromechanical  transducers  (gages) 
set  flush  into  the  surface  of  the  bodies.  The  experimental  body 
is  a 12- inch-diameter  hollow  dural  hemisphere  mounted  on  a suitable 
carriage  and  guided  during  a fall  of  11  feet  into  a tank  of  water 
by  two  vertical  rails.  The  velocity  at  impact  is  about  2k  fps. 

For  these  values  of  diameter  and  velocity  theifcV  pressure  to 
be  expected  would  be  about  1,600  psi,  and  have  a duration  of  about 
1/2  microsecond.  The  resonant  frequency  of  the  gages  used  here  is 
105  cps.  Thus  even  though  the  resonant  frequency  of  the  gages  is 
much  higher  than  that  of  the  pressure  plugs,  the  lower  entry  veloc- 
ities give  a pcV  pressure  duration  of  only  1/2  mi‘  rosecond.  For  a 
l/2-inch-diameter  gage  the  maximum  pressure  obtained  was  120  psi, 
as  compared  to  the  theoretical  value  of  1,600  psi.  However,  this 
experimental  value  of  pressure  was  registered  nearly  13  micro- 
seconds after  the  initial  impact,  which  means  that  the  peak  pres- 
sure might  have  been  modified  considerably  by  trapped  air.  It 
should  be  emphasized  that  the  experimental  work  was  done  with 
great  care  and  precision.  The  presence  of  trapped  air  seems  to 
be  indicated  by  this  experimental  procedure. 
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Appendix  A 

VON  KARMAN'S  DERIVATION 

A simple  derivation  for  the  impulsive  piston  pressure  p =(°cV 
is  given  by  Von  Karman  (Ref.  10).  When  a flat  plate  strikes  a water 
surface,  the  pressure  will  have  an  infinite  value  if  the  water  is 
considered  as  incompressible,  since  a finite  mass  of  water  is  given 
a certain  amount  of  kinetic  energy  in  z&xC-  time. 

It  is  possible  to  obtain  an  approximate  value  for  the  maximum 
pressure  taking  compressibility  into  account  in  the  following  man- 
ner: The  propagation  of  momentary  increase  of  pressure  in  a fluid 
takes  place  at  the  speed  of  sound  in  the  fluid,  designated  by  c. 
Therefore  the  mass  of  fluid  accelerated  in  the  time  dt  is<oScdt, 
where  S is  the  surface  of  the  fluid  struck  by  the  body.  Since  the 
velocity  of  this  mass  is  increased  from  zero  to  V in  the  time  dt, 
the  force  acting  is^ScV  and  the  pressure  is  p = (j’cV , 
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Appendix  B 

OWENS 1 TREATMENT  OP  THE 
TWO-SPACE  DIMENSION  PROBLEM 
BY  RIESZ'S  METHOD 


DEFINITIONS,  CONVENTIONS,  AND  THEOREMS 


The  Riesz  method  is  presented  in  detail  in  Ref.  6.  The  advan- 
tage of  this  method  is  that  it  shows  clearly  which  boundary  and 
Initial  conditions  are  needed  to  solve  a problem  in  hyperbolic 
equations.  As  a preliminary,  the  definitions,  conventions,  and 
theorems  necessary  for  the  application  of  Riesz* s method  arc  pre- 
sented. 

Definition:  Riemanr.-Liouville  Integral  (R-L  integral) 


l*f=  ty  = — [fCD/x-  ® [*■(*»)’>*''</* 

»'  rial  m Jo 

where  a is  a complex  number.  This  operator  associates  to  each 
function  f(x),  defined  for  x)>o,  a new  function  laf(x)  defined 
in  the  same  domain. 

Lemma : For  fixed,  bounded  f with  n continuous  derivatives, 
where  n Is  arbitrary,  and  for  fixed  x,  Iaf  Is  an  analytic  function 


of  a. 

Proof:  Iaf  as  defined  converges  for  He  ^ >0  and  is  analytic 
in  a.  So  when  Re  X.  integration  by  parts  yields 
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(Bl)  If  = 


f(o)  £ ^ / 

PI*  *■/')  Pfcti)  j0 


ff?1)(x-  fjVf 


However,  this  expression  equals  the  original  expression  in  Re  X'x) 
but  it  converges  for  /feh£>-/  * The  principle  of  analytic  contin- 
uation allows  the  use  of  Eq.  Bl  as  the  definition  of  Iaf  in 

> for  it  is  necessary  to  show  only  that  the  old  and  new 
expressions  are  analytic  in  a and  that  they  agree  in  a common  do- 
main- In  the  same  way  Iaf  can  be  defined  for  since  repeated 

integration  by  parts  will  allow  extension  of  the  domain  of  analyfc- 
icity  (1  unit  to  the  left  each  time)  throughout  the  complex  plane. 

The  use  of  analytic  continuation  is  the  basis  of  Riess's 
method  and  hence  the  reason  for  including  the  previous  proof. 

Properties  of  1°^: 


/,  Ilf~  f £(f)df 

* % A 

a,  If  - (■»-/)/ 7 / / 


g_ 

U* 


3.  I y — fU) 


4.  2 


< + J-\ 


J-! 


P J-  T iXJ 
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W 


f mi  d i 

J°  fy-p 


6.  l\  = X%*J 

if 


d*' 


8.  For  special  a,  i.e.,  a = 0 or  -n,  Iaf  depends  only  on  the 
values  of  f (f)  in  an  arbitrarily  small  neighborhood  of  f = x, 
but  for  general  a,  Iaf  depends  on  the  whole  range  of  values  of 


f ( f ) for  0 < f < X , 

Theorem : If  Hc6>Of  ne^i-  arbitrary,  then 
Application:  Abel's  integral  equation  is 


Ids'Ll*  ** 


fu)--  f%(V(x-i)~eJl  , o</i<! 

■Jo 

The  problem  is  to  find  g(x)  if  f(x)  is  known.  The  equation  can  be 
written 


,a  (t-ayi  /-£ 

fix):  p(i- /*.)-!—  f(t)(*-V  *ti  = nO-6)I  y 

PO-a)j° 


29 


NAVORD  REPORT  3323 


Since  l-$  >0  one  gets 

/ -0-&)  , 
— Z fa)  = -L 


no-s) 


ro-a) 


r--'A 

£ I f 


on  using  the  last  theorem  again  since  @>C>  • Prom  property  7 


go)  = 4-  r 


J_ ej__ 

PO-&)  <*X 


4 tL  f 


P(hA)d*  Pfe)J& 


Since 


r(/'4)  P(4)  = JEL 
s/t,  jr$ 


one  obtains 

3(x)  = 


5 - [km*-  t)  e'i/i 

7T  of  a Jo 


which  is  the  required  solution. 

Definition:  The  Riesz  integral  in  m independent  variables. 

c<~  tv, 

t t / £($)& 

J ' u /~j\  Jr, 


(B2)  f ) fiC  c*yv/}/<?){ 


fa)  Jp 


where 

f = f(x*,  x2,  . ..,  xm),  i.e.,  m independent  variables,  the 

superscripts  being  distinguishing  marks,  not  exponents. 

D = the  m-dlmensional  volume  bounded  by  the  m-1-diiuensional 
hypercone  with  apex  at  P(x^,  x2,  ...,  xm)  whose  equa- 
tion is  R = 0 and  x«—  % PC' 
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R = the  hyperbolic  distance  of  the  arbitrary  point 
tfC  f f **)  from  p?  given  by  R2  = ( f - /) 2 

-isx-xy-  — ( 

dVm  = d \ d "f < < o!  f "**,  an  m-dimensional  volume  element 
Hjjjt'a)  = a constant  defined  by 

!B3)  IU4  - e ) 


Remarks : If  depends  on  a,  f , P (corresponding  to  the  upper 
limit  of  uhe  R-L  integral)  and  on  a domain  D (corresponding  to  in- 
terval 0,  x of  the  R-L  Integral).  The  volume  P is  bounded  by  the 

t ? / 

hypercone  and  the  hyperplane  £ =0.  I$f  reduces  to  the  R-L  inte- 
gral for  m - 1.  Lastly,  the  Riesz  integral  is  m-fold,  but  by  in- 
troducing mean  values  for  f it  can  be  reduced  to  a double  integral. 
Theorem : 

(B4 ) s.  f(P)  = X^) 


In  particular  this  is  true  even  if  D is  a small  part  of  a cone  in- 
cluding the  apex  P.  This  is  analogous  to  the  R-L  integral  in  which 

I°f  depends  on  values  of  f ( f ) only,  in  a small  neighborhood  of 

1 = X / 

Theorem:  y 
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where  m(£,/|)  is  a mean  value  of  f.  The  essential  point  is  that 
this  formula  represents  I^f  in  terms  of  a simple  R-L  integral  over 
the  interval  (0,1). 


Theorem : 


where  D is  an  m-dimensional  volume  and  B is  an  (m  - 1) -dimensional 
volume,  i.e.,  the  "surface"  enclosing  D.  Differentiation  in  the 
co-normal  direction  is  represented  by  ^p,the  normal  being  the 
outward-drawn  normal  to  the  "surface."  The  m-dimensional  wave 
operator  defined  below  is  represented  by  O-yv),  This  is  Green's 
theorem  in  hyperbolic  geometry. 


If  the  components  of  the  normal  are  given  by 
7?  = ( A)  ^ij'"  > W-w)  the  co-normal  components  are  defined  by 
con  = ■"-r>  ■>*>  ) / This  follows  from  the  hyperbolic 

metric  associated  with  the  hyperbolic  differential  equation.  On 
the  plane  0,  the  normal  has  components  77  ) &) 

hence  con  = (~  f;  so  that  ^ = “ Jlf/  0n  the  Plane 

f 9)  ~Y)  =(&)~  ""J  provided  D (i.e.,  provided 

the  apex  of  the  cone  is  on  the  positive  side  of  the  plane  f ^=0  ) 
so  that  the  normal  extends  in  the  negative  direction  (see  sketch 
in  next  section).  Hence  con  = ( O j f ] and  ~ > 

The  expression  for  .sL  on  the  "curved  surface"  of  the  cone  will  not 
be  needed. 
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V) 


Recalling  that  in  m-dimensions  R^  = (% ~x')z—  j 

X + Z-yS)  Z J=Z  m -Z. 

one  can  compute  Oyy,  R-  E -2 — R.  — 2L  r~Te  ^ 

= o(,(«.-y/>l-Z)R  from  which  it  follows  that  for  a = 0,Q^ft  -zip 
and  hence  R^~m  satisfies  this  m-dimensional  wave  equation  (note 
that  Ora  is  defined  above).  An  attempt  to  build  up  a solution  from 
this  so-called  fundamental  solution  fails  when  m)>3  because  any  in- 
tegral involving  R2"m  diverges.  However  an  integral  of  Ra+2-m  where 
0C>w~4. will  converge,  and  its  limit  as  #->£>may  be  considered  as 
the  analytic  continuation  of  the  integral  of  Ra+2_m  into  the  domain 
^ yr)~4- } that  is,  it  defines  the  integral  of  R®  m.  This 
analytic  continuation  is  obtained  by  integration  by  parts  as  done 
in  the  continuation  of  the  simple  R-L  integral.  However,  this  in- 
volves a multiple  integral,  and  the  tool  available  is  Green's 
theorem,  the  use  of  which  corresponds  to  Integration  by  parts. 


TWC-SPACE  DIMENSIONAL  CASE 

Rlesz's  method  for  m = 3 will  be  applied  as  illustrated  in 
Pig.  12  to  give  physical  meaning  to  the  process.  The  case  m = 4 
will  be  treated  in  the  language  of  m = 3- 

Characteristic  cone,  vertex  at  P(ct,  y,  x): 

"Reflected"  cone,  vertex  at  P(ct,  - y,  x): 

cz(t-t)z- 

"Hyperbolic"  distances  from  P,  P to  Q(c\^)  fj  (Q  not  on 
cone)  are  respectively: 
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Let  Dp  Dp  be  the  respective  Interiors  of  the  two  cones  for  which 
V?  Let  s,  s be  the  surfaces  ^ -e>  and  interior  of  the 

cones ; s2  = Sg  the  surface  ty  = 0 and  interior  of  the  cones ; 
sp  s-l  the  remaining  "curved"  surfaces  of  tne  cones.  Let 
B~  S t S,/-S£  ) B-  5*  +5f  f-S'z  be  the  total  surfaces  of  the  two 
figures . 


34 


NAVORD  REPORT  3523 


The  two-space  dimensional  problem  is 
(37)  03^=  fii1r-CZ('PJi*+-Pyy)=4 

(B8)  X)  = O (initial  conditions) 

(B9)  %(t,0}  A)  = V}  UKX(-t)  5 <e(tjOjK)j  lxl>xte) 

(boundary  conditions) 

Now  Green's  theorem  is  applied  where  ^satisfies  U3  V?=  ^ and 

y=  r*-~'  ; r=  etu-0  ft*'3 

giving 

(bio)  «.(«.- -jjfh 

However 

yV=  — [?£*'& 

3 WJq 

from  Eq . BP  and 

/t3'ei  + Z)*  oC(eL-l)  HslX)  = z/rpfa-H ) 

from  Eq.  B3  and  from  the  properties  of  the  f'  -function  (duplication 
formula).  Thus  Eq.  BIO  becomes 

(bid  m -P(kn)r*v=Jj' <f\ h '§£)^s 

D 

Now  ±^if=  tffcjy)*)  from  Eq.  B4  so  that  Eq.  Bll  becomes 

..  / f /y^~- 
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Thus,  the  solution  f(t.  y,  x)  depends  only  on  the  values  of  ^ on 

the  boundaries  of  the  cone  D- . Since  R = 0 on  S-^,  and  with  o(>^) 

the  part  of  J on  S.  vanishes,  which  is  the  reason  for  not  express- 

ing  on  S,.  On  S,  — , - — — and  on  S_,  -2-  — 2-  . Then 
1 **  2 S>^ 

Eq.  B12  becomes 


<B»  > y'  *}  * & y*5hjf[~ 1 *** 


x-l 


t 


/ 


ffr, 


zrn(«n) 


4. 

73  £ 

i -X?-1  b%B,t 


l />\  rJfS 

L «*>■'»£ 


The  reflected  cone  will  now  be  used  but  the  solution  has  been 
written  in  the  form  of  Eq.  B13  to  call  attention  to  the  fact  that 
if  the  problem  had  no  boundary  conditions,  and  hence  the  cone  were 
not  cut  up  l j the  surface  ^=0,  the  Integral  over  S2  would  not 
occur  and  the  solution  ^ would  be  given  completely  by  the  initial 
conditions  which  appear  in  I / . In  this  case  the  reflected  cone 

JsJ 

would  not  be  needed.  However,  the  occurrence  of  mixed  boundary 
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conditions  renders  the  reflected  cone  a useful  device  enabling  one 
to  represent  the  solution  in  :erms  of  one  boundary  value. 

Since  the  initial  conditions  of  Eq.  B8  are  zero,  we  have 
from  Eq.  Bll 

(BiM 


Using  the  reflected  cone,  one  gets  the  analogous  relationship 

<b  if- 

3 JJ  & 7 


where 


= J. 

t—**  I »: 


f ^ a 

J4 


Prom  the  definitions  of  S2,  S2,  R,  and  R it  can  be  observed  that 
on  S2  = S2  where  # = 0,  R = R and  R^  = -R*.  Making  these  substi- 

i l 

tutions  in  Eq.  B15,  and  adding  Eq.  Bl4  and  Eq.  B15  gives 

zrnt**d(&  WWVMJ ? 


Lemma:  / t-0  that  is 




u* 


yy> 


/ 


oC^O  //-»/*€)  Jp 

Using  this  lemma,  Eq.  E17  becomes 


fit,  y,  <)  --  U*  ~ JJ*1 ' W 

$ 2 
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which  is  the  solution  to  the  problem.  For  a = 3 the  integral  con- 
verges so  that  by  carrying  out  this  liaising  process  directly 
yields 


Evaluation  of  this  integral  is  usually  very  difficult  and  is  not 
considered  here.  It  should  be  mentioned  that  in  this  evaluation, 
with  these  mixed  boundary  conditions,  the  method  developed  by 
J.  C.  Evvard  (Ref.  5)  may  be  used  when  the  point  P(ct,  0,  x)  is 
in  the  "subsonic"  region  and  only  when  y = 0. 


THREE-SPACE  DIMENSIONAL  CASE 


V 


Notation: 

V = entry  velocity  of  body 
c = velocity  of  sound  in  water 
p = density  of  water 
P = velocity  potential 
p - hydrodynamic  pressure 


Assumptions : 

1.  That  entry  velocity  remains  constant  and  — <<1. 

2.  That  c and  p remain  constant. 
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3-  That  flow  is  isentropic,  hence  inviscid. 

4.  That  effects  of  splash  and  of  gravitation  waves  are  ig- 
nored, and  the  water  surface  remains  plane. 

5.  That  the  body  is  very  blunt  at  the  point  of  impact:  hence 
the  boundary  conditions  are  applied  on  the  water  surface  and  not 
on  the  surface  of  the  body. 

6.  That  perturbation  velocities  of  the  fluid  are  small  so 
that  squares  of  these  velocities  are  neglected. 

The  problem  to  be  investigated  is  that  of  the  flow  pattern 
and  the  hydrodynamic  pressure  induced  by  the  body  entering  the 
water.  The  investigation  will  treat  only  the  first  few  moments 
of  contact  and  therefore  only  a small  penetration  of  the  body  into 
the  water.  Assumptions  4 and  5 are  made  for  this  purpose.  In 
particular,  the  pressure  history,  for  the  short  interval  considered, 
is  desired. 

Sine?  the  water  is  compressible,  tnls  treatment  differs  from 
the  classical  one  which  assumes  an  incompressible  fluid. 

The  flow  is  initially  irrotational  and  remains  irrotational, 
so  that  we  may  assume  a velocity  potential  f from  which  perturba- 
tion velocities  are  given  by  u = fx,  v = Py,  W = 

Prom  the  equations  of  hydrodynamics  which  are  linearised,  by 
neglecting  squares  of  velocities  and  using  the  assumptions,  the 
following  problem  is  obtained  whose  solution  answers  to  the  assump- 
tions and  description  »bov<?. 
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^ tt  " °2  (^xx  + ^yy  + ^zz)  =■  £7^  ^rC1 

p=  -Ke 


<P  (x,y,a,0)  = O' 
/t(x,y,s,0)  = 0. 


x2  + y2  + g2 


^ 0 (initial  conditions) 


Py(x,0,s,t)  = V j x2  + a2<  X(t)2} 

> (boundary  conditions) 
^(x,0,a,t)  =0  , x2  + z2  >X(t)2J 

The  expression  X(t)  describes  the  radius  of  the  ring  of  con- 
tact between  the  body  and  the  plane  water  surface  and  may  be  deter- 
mined from  the  geometry  of  the  body.  Hence  x2  + -a2  <C  X(t)-  corre- 
sponds to  points  on  the  body  in  the  disturbed  fluid  surface  and 
x2  + a2  >X(t)2  to  points  in  the  undisturbed  surface.  (Since 
boundary  conditions  are  applied  on  the  fluid  surface  rather  !••<? 
on  the  body  surface,  the  problem  is  that  of  an  expanding  disk  on 
the  surface,  sending  downward  pulses  corresponding  to  velocity  V.) 
The  following  definitions  are  needed: 

Characteristic  cone,  vertex  at  P(x,y,a,ct): 

xf- ( ?-y) z-  ~ o 


"Reflected"  cone,  vertex  at  P(x,-y,&,ct) : 

cHi -tf-  ( 1-xf -(?+?)*--  ti- if  = ^ 

"Hyperbol'c"  distances  from  P,  P to  Q(t>^f/v))  are  respec- 
tively  J- 

l c 1 i-rf -(</-? f-  (%-  2 

£ = f czCX-if-  tt-xf-  tt~  ^ 
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The  regions  Dj.,  Di,  S,  S,  S2  = S2,  Si,  "i,  B - S + Si  + 82  and 
B - + Si  + So  are  defined  as  in  the  two-dimensional  case  and 
are  described  later  (Fig.  14.). 

Applying  Green's  theorem  again  where 

and  . 

one  obtains 


<“»  « a*'**-/' 

D,  ->3 

but 


from  Eq.  B2  and 

#+«+*)-  « (°i-z)  H+fa ) = z**’ IT p(f  -h I ) P( t) 


from  Eq.  B3  and  from  the  properties  of  the /^-function.  Thus 
Eq.  B20  becomes 


L 


Now  take  a>3* 

- 0 , On  55,  2. 

5j/ 


Since  R = 0 on  Si  so  is  and  hence 

->  1 dJ> 

~-Si~  and  on  So,  - 2-  . Hence  Eq.  B2 
3C%  ~ 9 j? 


becomes 


HAVOFJD  REPORT  3523 


(B22) 


Since  the  initial  conditions  are  zero,  the  first  integral  vanishes 


This  process  applied  to  the  "reflected"  cone  gives 


/ 


(B2U) 


ffwUijl 


From  the  definitions  of  Sg,  S R*  R on  Sg  = Sg  where  j)  = 0,  R - R 
and  1^  -~^3.  Making  these  substitutions  in  Eq.  B24,  and  adding 


7 


Eq.  B23  and  Eq.  B24  gives 

(B25) 


- 2 


rVf  iV  = 


k si 


Using  the  lemma  that  I®  If  = 0 and  Eq.  B4,  Eq.  B25  becomes 
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(B2  6) 


_ / 


z«mi ^ 3 1 

where  the  "surface"  of  integration,  Sg;  is  the  volume  bounded  by 
the  surface  = 3 ' 

Again  it  can  be  seen  that  the  solution  is  completely  expressed 
in  terms  of  the  boundary  conditions  (and  initial  conditions  if 
different  from  zero).  In  this  case  the  integral  in  Eq.  B26  di- 
verges for  a = 0.  However,  the  factor  preceding  the  integral  was 
constructed  in  such  a manner  that  the  existence  of  the  limit  is 
insured. 

Evaluation  of  this  expression  (Eq.  B26)  for  is  quite  in- 
volved. The  special  point  P (0,0,0, ct)  in  the  supersonic  region 
can  be  treated  quite  simply  by  putting  R = 2(£~A)  ~ ”f  ““ J 2) 

introducing  polar  coordinates  in  the  (f;^)  plane,  using  the 
symmetry  of  Eif,  namely  ^ ( (?Co5>&j  £)j  fs/-n  &JX)=,f((,)‘)iJ  integrating 

with  respect  to  &■  , then  with  respect  to  £ and  finally  letting 
. However,  in  an  attempt  to  be  more  general,  mean  values 
will  be  used  to  express  ^(x,y,z,t)  for  the  general  point  P(x,y,z,cfc). 
For  this  purpose  put  H^foC-f-2)  = 2.  or  the 

factor  ^efore  the  integral  in  Eq.  B26  and  consider 


(n?7) 


_ ^ 
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Since  Sg  is  the  "surface"  ^ = 0,  introduction  of  the  $ -function 
makes  Eq.  B27  take  the  form 


(B28) 


[[-zftfA  S, 

>/ 


and  by  definition  from  Eq.  2,  this  can  be  written 


Now  Eq.  B5  provides  a formula  for  ijjf  with  which  one  can  discuss 
the  mean  value,  m((^  y ),  of  the  function  f mentioned  below. 

Referring  to  Pig.  14,  one  can  make  the  following  definitions: 

f~  ^ ( f ~~  X ‘ ) — Enc/i c/ ea  v d isfia-vee 

j-l 

from  / f j " > jf  Jto  the  "center  line"  of  the  cone  D in  the 

plane  "f  = constant. 

(B30) 


% = * 

\s~ 

-k 


V 


*z=  c*L  v) 2-  rz «.  3 ore i-  e*) 


sT  2 /\  , 


For  this,  problem  M = 4 and 

2_  / _ i/ 1 2.  , / i,  ,,\Z.  , / ■/?  -,\Z. 
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(B31) 


f 


\ ■=.  ct  -- 


and  Eq.  B2  becomes 


7 L0 

However,,/,  gives  the  function  inside  the  braces  of  Eq.  B32  with 
& / 

the  argument  replaced  by  one,  by  virtue  of  property  3 of  R-L  in- 
tegrals. Putting  (f  = 1 one  gets 


which  gives 
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If  the  cone  is  generalized  to  ra  dimensions , it  can  be  observed 
that  the  locus  of  points  r = constant, ^ = constant  represents 
the  (m~2) -dimensional  surface  of  an  (m-l)-dimensional  hypersphere 
con  iined  in  the  (ra-l) -dimensional  hyperplane  f = constant  whose 
center  is  at,  the  point  ( f X y"j  X 'Jwitk  /j  , 

Definition:  The  mean  value  of  f is  formed  by  integrating  f 
over  the  surface  of  this  sphere  and  dividing  by  the  area  of  the 
surface . 

By  induction  it  can  be  shown  that 


on-  / 


2?  2- 


represents  the  volume  of  an  (m-l) -dimensional  sphere.  This  for- 
mula holds  for  m = 3 (area  of  a circle)  and  m = k (volume  of  an 
ordinary  sphere.) 

Now  assume  that  the  formula  is  valid  for  an  (m-l)-dimensional 
sphere  and  prove  that  it  holds  for  an  m-dimensional  sphere.  By 
induction  it  will  then  hold  for  all  dimensions.  By  the  volume  of 
an  m-dimensional  sphere  is  meant  the  "area"  of  any  cross  section 
multiplied  by  the  elementary  thickness  and  Integrated  throughout 
this  sphere  (see  Fig.  15).  The  area  of  this  cross  section  is  the 

volume  of  the  (m-l)-dimensionnl  sphere.  Thus 

/7[  , Tr>-/ 

* ULJ  (rc*se)c/& 

jo  •»>-/ 
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Using 

(B34) 


which 

sphere 

B, 

ume  (ai 
the  (m 

(335) 


(B3*0 


which  has  the  same  form  as  the  formula  for  an  (m-l)-dimensional 
sphere  and  the  assertion  is  proven. 

By  considering  the  area  as  the  radial  derivative  of  the  vol- 
ume (and  putting  m-1  for  a in  Sq.  B34  one  gets  for  the  area  of 
the  (m-2)-dlmensional  surface  of  an  (m-1 ) -dimensional  sphere 


(B35) 


'yyj-Z  P/T  z 
^ /~7  / -^  - / j 
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and  by  the  above  definition 


(B36)  = 


/ 

5yyi-cL 


3 fixed 


where  dw  is  the  surface  element  of  the  (m-2) -dimensional  surface 
of  the  sphere  whose  center  is  the  point  (A"7— A;Xy'/  j Xnn)  • For 
the  present  case  W = p ~ 1 (implying  r = .A  from  Eq.  B30)  and 
the  center  is  the  point  (x}y}tjCt-  % ^),Then  Eq.  B36  becomes 


W//-  2)- 


j--  j TLfyjdw  =4yFz Jffl  fj  %) 

/"-A  $!xed 


where  use  has  been  made  of  Eq.  B35.  Substituting  this  in  Eq.  B34 
gives 


Cir  _ 

(B37)  //=  j-  f 42.  f 

* 4/rJ  J J_ 

o 


Applying  this  formula  to  Eq . B29  where  -f  = ~ 2 
the  solution  given  by  Eq.  B24  then  becomes 


SjfY '?,^?,i)SC0 
3J 


Cir 


(B38) 


^ \ sUt(f,o>%,ciz5)S0i)tlvi/ 

Tjjy  * 

r*t 
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Appendix  C 

INSTANTANEOUS  PULSE  STRENGTH 


The  potential  equation  is 


Following  Lagerstrom,  pc  i coordinates  are  introduced  (Fig.  l6) 


and 


hi 

x-f 
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For  fixed  z 

r,  Jr,  = e^e 


hence  the  element  of  area  is 


c/4  = - r,  dr, 


The  potential  equation  becomes 

Ct  ZF 


LP(. ,y, J J 


2 D 

Differentiating  with  respect  to  2r 

ZF  , Cir 


-aj£(x,y,i,t)  - j Je  2-  i f(e)91i-Z-)c/> 7 


Now,  from 


r'se'+i*  , £f  =-! 

•l  " ' 11  I (7 

zr  , ct 


1 
</ 


- PPi. 


d 2- 


n 


Therefore 


«=//  . -t-  -j  v 


but 


P 2-  3 2- 


if 

32- 

?£ 

<? 

5(? 

Hence 


j:T  ct  ?F 

* * 


im 


Setting  ?:  = 0 one  obtains 
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y*  (*,y,  o,  t)  = - ztt  o,  t) 


Tine  Instantaneous  pulse  strength  is  thus  directly  proportional 
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Appendix  D 

LEGENDRE  STANDARD  FORM 


The  Legendre  standard  form  of  the  elliptic  integral  of  the 
first  kind 

If 

r f de 

(di  ) l « / 7==  ="■  - ■■ 

/?]/( //-  r'cose)  z+-(t‘ 


This  equation  can  be  reduced  to  the  Legendre  standard  form  by  the 
following  four  transformations: 

Transformation  1 

Let  % ~ CD 5,  «$> 

Upon  substitution  into  Eq.  Dl  and  factorization  of  the  expression 
under  the  radical  one  obtains 


(D2) 


i- i T)  ’ 


-lfJL 

r'Jnw 


where 

* = - jrj/  v-  L JiCpr) 
g=-jr,(l-  * 

Y=  / , S = -/ 
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Transformation  2 

To  rid  the  radical  of  odd  powers  of  the  variable  of  Integra- 

/ / . * 

tion,  the  procedure  outlined  by  Karman  and  Biot  (Ref.  11/,  for  a 
fourth-degree  polynomial  is  followed,  and  the  homographic  trans- 
formation made 


% - t 

/ 1-  * 

This  yields  the  integral 

r = -L  / (?- 

r'  1 — — — . 

J V-  0,(2) 

where  -x)/f  -/$)lzJ[(p-?yp-  0+ 

and  the  requirement  that  odd  powers  of  vanish  determines  f and 

q as  follows 

f « - cl  - V*2-/ 

$■  ~ ~ + V 


where 

a = L±l! 

2r' 

Upon  factoring  one  can  rewrite  the  integral  as  follows 


(D3)  £ = K 


f 

J i-Tit-  5?  ‘i  l)(ri~f(  i y ' 


where 


K = 


r'-fcp- *.)(?-  a)  ( f-  *)(p- 1)  ‘ 
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3 


cl 


(P-  v)  { p~  <) 

Li±; 

a rt 


-y^-77  i—  ci. 
-jcf-j1  — cl 

-j^P  - g 

-m1  f-  % 


Transformation  3 

Since  a.'y  I and  cL>C>  ir ^ /J) therefore  If-2/  > I h2) , 

Following  Karman  and  Blot,  putting 


where  C<  j ; hence,  Eq.  D3  becomes 
(D4)  — r 

3 J i/l^nxrrW) 

Transformation  4 

Lastly,  applying  the  transformation 


f - 


/ 


Equation  D4  becomes 

K. 


4/ ^ 


S-fPTc^  J -)/- ^ s;r>cf 


or,  with  the  limits  of  integration 
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(D5)  I = 


*&',**)  [ 

Jo 


f 

-)//  - A.  t ')  s/v2  r i 


where 


2 

-tar')* 7 £ 


H-  t'~  zr'  ^ 


Equation  D5  is  a complete  elliptic  integral  of  the  first  kind. 
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NOMENCLATURE 


A Area 

a Radius  of  sphere 
C,  c Velocity  of  sound  in  water 

P Surface  of  fluid  struck  by  the  body 
F,  F'  Points  on  the  -nit  circle  of  Pig.  2 
f Instantaneous  pulse  strength 
m Mach  master 
m Number  of  dimensions 
P Field  point 

?,  U,  V Functions  of  the  complex  variable 
P,  u,  v 

p Pressure 

R Hyperbolic  distance 

R(t)  Instantaneous  radius  of  the  disk 

R(t)  Rate  of  growth  of  radius  of  wetted  area 

r A point  on  the  disk 

S Wetted  surface  of  the  body 

S'  Horisontal  free  surface  of  the  liquid 

t Fixed  time 

tc  Time  at  which  rate  of  expansion  of  wetted  surfaces  Is 
sonic 
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Velocity  vec^r  in  flow  field 
Normal  velocity  of  disk 

„ hJe. 

/ 

Coordinate 

Entry  angle 

J7 

Coordinate  axis  of  the  plane 
7,  h,  ^ Coordinates  of  a variable  point  in  space 
<p  Density  of  water 
O' , X Polar  coordinates 
X Variable  time 
f Velocity  potential 

■y  Angle  between  HD  and  BF  in  Fig.  2 


U 

V 


x 

x#  f 


9 

A 


f/-? 
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